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ABSTRACT 

Merger-remnant galaxies with kiloparsec (kpc) scale separation dual active galactic nuclei (AGNs) 
should be widespread as a consequence of galaxy mergers and triggered gas accretion onto supermassive 
black holes, yet very few dual AGNs have been observed. Galaxies with double-peaked narrow AGN 
emission lines in the Sloan Digital Sky Survey are plausible dual AGN candidates, but their double- 
peaked profiles could also be the result of gas kinematics or AGN-driven outflows and jets on small 
or large scales. To help distinguish between these scenarios, we have obtained spatial profiles of the 
AGN emission via follow-up long-slit spectroscopy of 81 double-peaked narrow-line AGNs in SDSS 
at 0.03 < z < 0.36 using Lick, Palomar, and MMT Observatories. We find that all 81 systems 
exhibit double AGN emission components with ~kpc projected spatial separations on the sky (0.2 
h^Q kpc < Ax < 5.5 hjQ kpc; median Ax = 1.1 hj^ kpc), which suggests that they are produced 
by kpc-scale dual AGNs or kpc-scale outflows, jets, or rotating gaseous disks. Further, the objects 
split into two subpopulations based on the spatial extent of the double emission components and the 
correlation between projected spatial separations and line-of-sight velocity separations. These results 
suggest that the subsample (58^ %) of the objects with spatially compact emission components may 
be preferentially produced by dual AGNs, while the subsample (42^ %) with spatially extended 
emission components may be preferentially produced by AGN outflows. We also find that for 32^g % 
of the sample the two AGN emission components are preferentially aligned with the host galaxy 
major axis, as expected for dual AGNs orbiting in the host galaxy potential. Our results both narrow 
the list of possible physical mechanisms producing the double AGN components, and suggest several 
observational criteria for selecting the most promising dual AGN candidates from the full sample of 
double-peaked narrow-line AGNs. Using these criteria, we determine the 17 most compelling dual 
AGN candidates in our sample. 

Subject headings: galaxies: active - galaxies: interactions - galaxies: nuclei 



1. INTRODUCTION 

It is now established that most bulge-dominated galax- 
ies host central supermassive black holes (SMBHs) and 
that galaxies frequently merge with one another. As a 
natural consequence, there must exist a subpopulation 
of galaxies that harbor two SMBHs. These SMBH pairs 
are observable if they are accreting gas as active galactic 
nuclei (AGNs), though some may be enshrouded by dust. 
The kiloparsec (kpc) scale separation dual AGNs are 
an intermediate evolutionary stage between the galaxy 
pairs separated by tens of kpc use d as proxies for the 
merger rate (e.g., prLejal] [20081; lEUison et all 120081 : 
Ide Ravel et alll201lh and the sub-pc separation binary 
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SMBHs (e.g., Er acleous et al.l feoil) that are expected to 
produce gravitational waves upon coalescence. Although 
there are thousands of known galaxy pairs with tens of 
kpc separations , only a handful of kpc-scale dual AGNs 
(iKom ossa et al. 2003; Hudson et al. 2006t iBianchi et al.l 
120081: iKoss et al.ll201lb iFu et aLll2011bh , one 100 pc scale 
AGN pair dFabbiano et al.l [201 lh . and one 10 pc scale 
AGN pair ([Rodriguez et al.l 12006) have been confirmed. 

The small number of confirmed dual AGNs is in con- 
flict with expectations based on the galaxy merger rate 
and the triggering of gas accretion o nto SM BHs dur- 
ing mergers (e.g., |Sp ringel et al. 120051 : IHopkins et al.l 
120051 : IVan Wassenhove et al.l l2012f ). Although discov- 
eries of individual dua l AGN systems co ntinue (e.g., 
iComerford et al. 2009ri IB arrows et al.ll2012h . a system- 
atic survey of dual AGNs is now necessary to close the 
gap between the observed and expected numbers. 

Such a systematic survey of dual AGNs would shed 
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light on drivers of galaxy evolution such as galaxy merg- 
ers and SMBH growth. As direct observational tracers 
of galaxy mergers, dual AGNs can be used to estimate 
the galaxy merger rate. Further, the number of dual 
SMBHs powering AGNs can set constraints on the acti- 
vation probability of SMBHs in merging galaxies, which 
is a key parameter in models of the growth of SMBHs by 
gas accretion during galaxy mergers. 

With the advent of large spectroscopic surveys of 
galaxies, it is now possible to conduct a systematic cen- 
sus of dual AGNs. Two dual AGN candidates were 
identified in the DEEP2 Galaxy Redshift Survey based 
on their long-slit spectra, which revealed two compo- 
nents of AGN-fueled [O III] A5007 emission separated by 
~kpc spatially on th e sky and a few hundred km s -1 in 
line-of -sight velocity ( Ger ke et al.ll2007t IComerford et al.l 
2009a). The velocity separations manifest as double- 
peaked emission lines in the one-dimensional spectra. 
However, the extreme rarity of double-peaked AGNs de- 
mands a larger survey for statistical studies. 

Recently, 340 unique AGNs with double-peaked 
[O III] A5007 were identified in the Sloan Digital Sky 
Survey (SPSS] at 0.01 < z < 7 (|Wang et al.l [20091 : 
iLiu et al.ll2010bl : ISmith et al.ll2010f ). These are plausi- 
ble dual AGN candidates, but velocity offsets in AGN 
emission lines can also be produced by narrow-line region 
(NLR) effects such as small-scale (< 100 pc) gas kinemat- 
ics an d outflows (e.g. JDas et al.ll2005l l2006t ISmith et al 
2012) or larger-scale AGN outflows (e.g. JCrenshaw et al 



201ft iFischer et al] l20TTh . Since SDSS fiber spectra 



carry insufficient spatial information, spatially resolved 
follow-up slit spectroscopy is required to determine both 
whether the double-peaked emission lines come from spa- 
tially distinct regions and what the physical scales of the 
regions are. These observations will help constrain which 
physical mechanisms are producing the double-peaked 
AGN lines and which objects are good candidates for 
dual AGNs. 

We have obtained follow-up long-slit observations of 81 
(or one-fourth) of the double-peaked narrow-line AGNs 
in SDSS using a combination of the Kast Spectrograph 
at Lick Observatory, the Double Spectrograph at Palo- 
mar Observatory, and the Blue Channel Spectrograph at 
MMT Observatory. Since the orientation of the double 
emission components on the sky is unknown, we observe 
each object at two position angles, generally separated 
by 90°, so that we can determine the full spatial separa- 
tion of the two emission components on the sky as well 
as their position angle. We use this information, as well 
as the spatial extent of the emission and existing multi- 
wavelength observations of the host galaxies, to constrain 
what mechanisms (gas rotation, AGN outflows and jets, 
or dual AGNs) on what scales (< 100 pc, ~kpc, or ~10 
kpc) produce the line profiles, and to identify which ob- 
jects are the best candidates for dual AGNs. 

We assume a Hubble constant Ho = 70 km s _1 Mpc -1 , 
Vtm = 0.3, and = 0.7 throughout, and all distances 
are given in physical (not comoving) units. 

2. THE SAMPLE AND OBSERVATIONS 



2.1. The Sample 

We selected targets from the three catalogs that 
have identified active galaxies in SPSS with double - 
peaked [O III] A50 07 emi ssion lines: I Wang et al.l (j2QQ9l ). 
ILiu et all (|2Q10b[ ). and ISmith et al.l (l2Q10ft (we note 
that one of the o b jects in ISmith et al.l 120101 was origi- 
nally identified in |Xu^Eomossil2009|) . These AGNs 
were selected usin g; sta ndard emission line diagnos- 
tics ([Baldwin et al.l 119811 : iKewlev et aTll200lD . and the 
double- peaked AGN were selected with varying tech- 
niques. Wang et "aD (j2QQ9h select 87 double-peaked Type 
2 AGNs that have blueshifted and redshifted compo- 
nents of [O III] A5007, relative to the host galaxy red- 
shifts, that have a wavelength difference > 1 A and 
a flux ra t io betw een the two components of 0.3 - 3; 
ILiu et all (|2Q10bf ) select 167 double-peaked Type 2 AGNs 
by eye, through visual identification of systems that 
have similar doubl e-peaked pro files in both [O III] A5007 
and [O III] A4959; ISmith et al.l (|2010f ) select 86 double- 
peaked Type 1 AGNs and 62 double-peaked Type 2 
AGNs through visual selection of systems with double 
peaks in both [O III] A5007 and [O III] A4959. After ac- 
counting for the overlaps between these samples, there 
are 340 unique double-peaked AGNs. We drew our tar- 
gets from these 340 objects. 

Since our aim was to use follow-up long-slit spec- 
troscopy of these objects to measure ~kpc spatial off- 
sets between the [O III] A5007 emission components, we 
selected targets as double-peaked AGNs that had 1) air- 
masses below 1.5 during our observing runs, 2) redshift s 
z < 0.36 where <kpc spatial offsets were resolvable with 
the spectrograph pixel scale, and 3) sufficiently bright 
r-band magnitudes r < 18.5 to minimize exposure times 
and maximize the number of objects observed in this 
campaign. The latter two criteria are spectrograph- and 
telescope-dependent, and we selected targets for each 
telescope based on the spectrograph pixel size and the 
aperture of the telescope as described below. We made 
no selections on any other galaxy properties. 

We obtained follow-up slit spectroscopy for a sample 
of 81 double-peaked AGNs at 0.03 < z < 0.36 and 
14.5 < r < 18.5. These targets are a representative 
sample of the 237 double-peaked AGN in SDSS that fit 
the criteria z < 0.36 and r < 18.5; Kolmogorov-Smirnov 
tests find a 91% (87%) probability that the redshifts (r- 
band magnitudes) of our observed sample and the 237 
potential targets were derived from the same parent dis- 
tribution. We note that this is a pilot program, and long- 
slit observations of more double-peaked AGN in SDSS 
are forthcoming (J. M. Comerford et al., in preparation). 

2.2. Observations 

We obtained long-slit spectroscopy of the 81 double- 
peaked AGNs using the Kast Spectrograph at the Lick 3 
m telescope (pixel size 0'/78), the Double Spectrograph 
(DBSP) on the Palomar 5 m telescope (pixel size 0'/47 
for the red detector and 0'/39 for the blue detector), and 
the Blue Channel Spectrograph on the MMT 6.5 m tele- 
scope (pixel size Of! 29). We optimized the aperture size 
and pixel size of each telescope-spectrograph combina- 
tion to observe the lower redshift and brighter objects 
with Lick/Kast, the medium redshift and brightness ob- 
jects with Palomar/DBSP, and the higher redshift and 
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TABLE 1 

Summary of Observations 



SDSS Designation 


Telescope /Instrument 


Observation Date (UT) 


01 (°) 


02 (°) 


Exposure Time (s) 


SDSS J000249.07+004504.8 


Lick/Kast 


2009 August 17 / 2009 August 18 


69.4 


158.0 


3600 


SDSS J000656.85+154847.9 


MMT/Blue Channel 


2010 November 6 


44.4 


134.4 


1080 


SDSS J000911. 58-003654.7 


MMT/Blue Channel 


2010 November 6 


67.0 


157.0 


1080 


SDSS J011659.59-102539.1 


MMT/Blue Channel 


2010 November 5 


28.8 


118.8 


1080 



Note. — We observed each object at two position angles, 6\ and 62 (given in degrees east of north), and the exposure time 
given is for each position angle. 

(This table is available in its entirety in a machine-readable form in the online journal. A portion is shown here for guidance 
regarding its form and content.) 



TABLE 2 

Summary of Measurements 



SDSS 


Redshift 


Spec. 


Av 


Ax 


Ax 


Phys. 


#sky 


AO 


e 


Other Obs. 


Reference 


Name 




Type 


(km s- 1 ) 


(") 


(h~Q kpc) 


Extent 


(°) 


(°) 








0002+0045 


0.087 


2 


511 ± 3 


0.58 ± 0.05 


0.95 ± 0.08 


c 


63.9 ± 5.0 


46.3 


0.106 


NIR/SS/VLBA 


1, 2, 3, 4, 5 


0006+1548 


0.125 


2 


359 ± 3 


0.37 ± 0.02 


0.84 ± 0.05 


e 


172.6 ± 3.9 


51.8 


0.040 


NIR 


2, 3 


0009-0036 


0.073 


2 


304 ± 5 


0.25 ± 0.01 


0.35 ± 0.01 


c 


56.2 ± 1.1 


79.1 


0.188 


NIR/SS 


1,4 


0116-1025 


0.150 


2 


288 ± 5 


1.02 ± 0.06 


2.68 ± 0.16 


c 


117.0 ± 2.6 


1.8 


0.151 


NIR/SS 


1, 4 



Note. — Column 3 shows the spectral type of the AGNs. Column 4 shows the velocity offset between the [O III] A5007 peaks in the SDSS 
spectrum. Column 5 shows the angular projected spatial offset between the two [O III] A5007 emission features on the sky, as measured from 
our slit spectroscopy. Column 6 shows the physical projected spatial offset. Column 7 shows whether the double AGN emission features in 
the two-dimensional spectrum appear compact (c) or extended (e) in physical extent. Column 8 shows the position angle between the two 
[O III] A5007 emission features on the sky, measured from our slit spectroscopy in degrees east of north. Column 9 shows the absolute value 
of the difference between the measured position angle and the isophotal position angle of the major axis of the object from SDSS r-band 
photometry (with the same error as in Column 8). Column 10 shows the ellipticity of the object from SDSS photometry. Column 11 lists other 
published observations of the object: near infrared imaging (NIR), slit spectroscopy (SS), Very Long Baseline Array observations (VLBA), 
integral field unit spectrosc opy (IFU), Hubble Spac e Telescope imag ing (HST ), Very Large Arra y observations (VL A), and Chandra observation s 
(Chandra). References: (1) Liu et al. (2010b), (2) Wang et al. (2009), (3) Fu et al. (2012), (4) Shen et al. (2011), (5) Tingav & Wavth (2011), 
(6) Smith et al. (2010), (7) Rosario et al. (2011), (8) Fu et al. (2011a), (9) McGurk et al. (2011), (10) Rosario et al. (2010), (11) Comerford et al. 
(120111 ). 

(This table is available in its entirety in a machine-readable form in the online journal. A portion is shown here for guidance regarding its form 
and content.) 



fainter objects with MMT/Blue Channel. We note that 
since we resolve each double-peaked AGN's two emission 
components in velocity space, it is possible to measure 
their positional centroids to much better than the seeing 
limit and that the limiting angular scale is instead the 
pixel scale. 

We observed subsamples of our target list as follows: 
16 objects at 0.04 < z < 0.16 (median z = 0.08, where 
1" corresponds to 1.5 h^Q kpc) and 14.5 < r < 17.5 
with Lick/Kast; 16 objects at 0.06 < z < 0.15 (me- 
dian z = 0.11, where 1" corresponds to 2.0 h^Q kpc) 
and 15.0 < r < 18.2 with Palomar/DBSP; and 49 ob- 
jects at 0.03 < z < 0.36 (median z = 0.15, where \" 
corresponds to 2.6 hj^ kpc) and 14.5 < r < 18.5 with 
MMT/Blue Channel. Table [1] summarizes the observa- 
tions. We used a 1200 lines mm -1 grating with each tele- 
scope/spectrograph, centered such that the wavelength 
range spanned H/3 and [O III] for the redshift range of 
each night's subsample of targets. 

We observed each target twice, with the slit at two 
different position angles, in order to determine the orien- 
tation of the AGN emission components on the plane of 
the sky. We typically made observations at the isophotal 
position angle of the major axis of the object in SDSS r- 
band photometry and the corresponding orthogonal po- 
sition angle. However, if the SDSS imaging showed a 
companion near our target we aligned the position angle 



to include that companion rather than using the position 
angle of the galaxy, and at times mechanical constraints 
on the rotation of the telescope prevented us from ob- 
serving at a second position angle that was orthogonal 
to the first. 

The data were reduced following standard procedures 
in IRAF and IDL. 

3. ANALYSIS OF SDSS DATA 
3.1. Velocity Measurements 

In this section we present measurements of the veloc- 
ity offsets between the redshifted and blueshifted compo- 
nents of the double-peaked [O III] A5007 emission lines 
in the SDSS spectra. Althoug h velocity offsets have been 
measured in the literatur e ([Wang et al.l 120091 : iLiu et all 
I2010bl : ISmith et al.l 12010). cases where different authors 
measured the velocity offset for the same object often 
conflict. Consequently, for completeness and consistency, 
we measure all of the velocity offsets from the SDSS spec- 
tra here (Table EJ). 

For each galaxy, we fit two Gaussians to the 
continuum-subtracted [O III] A5007 emission line profile, 
which is the emission line with the highest signal-to-noise 
ratio in our observations. We measure the line-of-sight 
velocity difference based on the wavelengths of the peaks 
of the best-fit Guassians, and the error in the velocity 
difference is derived from the errors in the peak wave- 



4 



Comerford et al. 



lengths of the best-fit Guassians added in quadrature. 

3.2. Ellipticities 

Using the r-band Stokes para meters Q r and U T from 
the SDSS photometric pipeline ([Stoughton et al. 2002), 



we compute the ellipticity e = yjQ^-^U^ for each object 
in our sample. While adaptive moments can be used for 
more refined measurem ents of galaxy shapes and are used 
in weak lensing studies ([Bernstein fc Jarvisll2002[ ) , for our 
purpose of measuring the overall ellipticity of a galaxy's 
light the Stokes parameters are the most appropriate. 

4. LONG-SLIT SPECTRAL ANALYSIS 
4.1. Spatial Separation Measurements 

For each two-dimensional long-slit spectrum, we deter- 
mine the projected spatial separation between the two 
[O III] A5007 emission features by measuring the spatial 
centroid of each emission component individually. We 
measure the spatial centroid of an emission component 
by first centering a 2 A (rest-frame) wide window on the 
emission. Next, at each spatial position we sum the flux, 
weighted by the inverse variance, over all wavelengths 
within the window. We place a 10 pixel window around 
the emission, and define the window's center to be the 
spatial position where the summed flux is maximum. We 
then fit a quadratic to the summed flux to locate a peak, 
define a narrow window centered on the peak flux, and 
compute the line centroid within this window. We derive 
the error on this spatial centroid by repeatedly adding 
noise to the spectrum drawn from a Gaussian with the 
variance of the pixels in the window and redoing all cen- 
troid measurements. 

For a given target, we measure the projected spatial 
separation between the two [O III] A5007 emission com- 
ponents at both position angles observed, and then com- 
bine the two spatial separations to yield the spatial offset 
between the two emission components on the sky and the 
position angle of their offset. If observations at position 
angles 9\ and 62 yield spatial separation measurements of 
x\ and X2, respectively, then the true position angle # s k y 
of the two emission components on the sky is determined 
through the numerical solution of 

xi cos(# sky - 2 ) = x 2 cos(<9 sky - 0i ) . (1) 

With # s ky determined above, the spatial separation of 
the two emission components on the sky is then Ax = 
x\j cos(# s ky — #1) or equivalently, Ax = £2/ cos(# s k y — 62)- 

4.2. Spatial Extent of the Emission 

We also classify by eye whether the AGN emission com- 
ponents in the long-slit spectra are spatially "compact" 
or "extended". We find that each emission feature dis- 
plays two components, and if both components appear 
spatially distinct and compact at both position angles 
observed we label the object "compact" . If one or both 
of the components appear spatially extended or diffuse 
at one or both position angles, we label the object "ex- 
tended". This type of extended structure is known to 
be produced by gas kinematics, outflo ws, or jets (e.g., 
iRosario et al.ll2010[ iFischer et al.ll20llh . though we note 
it could also be produced by dual AGNs where one or 
both of the AGNs also has visible emission from gas kine- 
matics, outflows, or jets. 
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Fig. 1. — Segments of the two-dimensional long-slit spectra 
and one-dimensional SDSS spectra of five example double-peaked 
AGNs that exhibit two spatially compact emission components. 
Each object was observed at two roughly orthogonal position an- 
gles, as labeled in the two-dimensional spectra. The spectra are 
shifted to the rest frame of the host galaxy, and the dotted ver- 
tical lines denote the rest wavelengths of H/3, [O III] A4959, and 
[O III] A5007. Night-sky emission features have been subtracted 
from the two-dimensional spectra. 
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Fig. 2. — Same as Figure [TJ but for three example double-peaked 
AGNs that exhibit spatially extended emission. 

We find 47 spatially compact objects and 34 spatially 
extended objects, corresponding to 58^6 % an d 42^ % of 
the sample, respectively. Figure Q] shows examples of the 
spatially compact object spectra, while Figure [2] shows 
examples of the spatially extended object spectra. 

Sections [5] and [6] address other observations of our sam- 
ple and some objects that we find have unusual long-slit 
spectra; for results for the complete sample skip to Sec- 
tion H 

5. MULTIWAVELENGTH OBSERVATIONS FROM 
THE LITERATURE 



In an effort to classify the nature of the double-peaked 
[O III] A5007 emission lines in SDSS spectra, some au- 
thors have conducted follow-up observations of subsam- 
ples of these objects in optical, near-infrared (NIR), ra- 
dio, and X-ray. These observations aimed to test the dual 
AGN hypothesis by resolving systems with two AGNs 
or two stellar components, where the stellar components 
are presumed to be remnants of the progenitor galaxies 
from the merger that produced the dual AGNs. Here we 
summarize these observations and their overlap with our 
sample (see also Table [2]). 

The follow-up optical slit spectroscopy and NIR 
imaging of 31 Type 2 Seyfer t s with double-pe aked 
[O III] A5007 (jLiu et al.l l2010at IShen et al.l 12011ft in- 
cludes 14 objects in our sample. None of our objects 
exhibit two stellar components in the NIR imaging. How- 
ever, in three of our objects we measure that the emis- 
sion components in the long-slit spectra are separated by 
less than the ~ 0'/4 spatial resolution of the NIR imag- 
ing, such that stellar nuclei coincident with the emission 
components would not be spatially resolved in the NIR 
imaging. In the remaining 11 objects, stellar nuclei that 
are spatially coincident with the double emission com- 
ponents could be resolved in the NIR imaging, so either 
they have no secondary stellar component at these spatial 
separations, or the secondary stellar nucleus is obscured 
or too faint to be visible in the NIR imaging. 

Follow-up adaptive optics (AO) imaging strived to re- 
solve smaller-separation stellar compon ents in 55 galax- 
ies with double-peaked [O HI] A50 07 (|Fu et alJ l2011at 
iRosario et all 120111 : iFu et al.ll2012h . In addition to ob- 
taining integral field unit (IFU) spectros copy of 42 galax- 
ies with double-peaked [O III] A5007, IFu et al.l (|2012h 
also analyzed high-resolution Advanced Camera for Sur- 
veys/ WFC2 archival Hubble Space Telescope (HST) im- 
ages of six additional double-peaked galaxies. These 
IFU, AO, and HST observations include 34 objects in 
our sample, of which seven show double stellar compo- 
nents in AO. 

Although position angles between th e stellar compo- 
nents are not reported in the litera ture (|Fu et al.ll2011at 
IRosario et al.ll201 J ; IFu et all 12012), we estimate the po- 
sition angles for these seven objects based on our best 
estimates of the bright centroid of each component in 
the published images. In Section l7.4.2[ we compare these 
position angle estimates, as well as the reported spatial 
separations, to the position angles and spatial separa- 
tions measured from our long-slit spectroscopy. 

Further, Very Long Baseline Array (VLBA) observa- 
tions of 11 Type 2 Seyferts that have double-peaked 
[O III] A5007 as well as Faint Images of the Radio Sky 
at Twenty-cm (FIRST) detecti ons reveal compact radi o 
emission in only two objects (Tingav & Wavthl [201 lft . 
One (SDSS J151709.21+335324.7) of these two galaxies 
is part of our sample, and this galaxy also has a radio 
jet that suggests a radio jet-driven outflow produces th e 
double-peaked [O III] A5007 lines ([Rosario et all 12010ft . 
None of the VLBA observations show double radio cores, 
although Expanded Very Large Array observations of an- 
other double-peaked AGN (not in our sample) revealed 
two radio cores t hat confirmed it as a dual AGN system 
(|Fu et al.ll20lTbh . We also find no significant correla- 
tion between FIRST detections and other properties of 
double-peaked AGNs that are measured here, such as 
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Fig. 3. — Same as Figure [TJ but for three noteworthy galax- 
ies in our sample that have unusual emission features. SDSS 
J080740. 99+390015. 2 (top) has triple-peaked narrow lines and 
single-peaked broad lines; SDSS J000911. 58-003654.7 (middle) has 
a companion spiral galaxy that may be at a similar redshift; and 
SDSS J100145. 30+283330. 3 (bottom) has spatially extended emis- 
sion linking it to its companion galaxy, which is a possible signature 
of gas outflow or inflow. 



line-of-sight velocity separation and projected physical 



Fig. 4. — SDSS images of the three galaxies whose spectra 
are shown in Figure SDSS J080740. 99+390015. 2 (top) may 
host AGN outflows or jets; SDSS J000911. 58-003654.7 (middle) 
has a companion spiral galaxy visible to the southwest; and SDSS 
J100145. 30+283330. 3 (bottom) has a companion galaxy visible to 
the northwest. 

separation. 

Finally, Chandra observations of the double-peaked 
AGN in SDSS J171544.05+600835.7 show two X-ray 
sources with separation Ax = 1.9 hj^ kpc and 
orientation # s k y = 147°, which are consistent with 
the separation and orientation measured between the 
double [O HI] A5QQ7 s ources in the long-slit spectra 
(|Comerford et al.ll201l[ ). While this is strong evidence 
for dual AGNs, additional follow-up observations are re- 
quired for confirmation of dual AGNs in this system. 

6. UNUSUAL SPECTRA 

Three of our targets have particularly unusual and 
noteworthy spectra, which we discuss briefly here. Fig- 
ure [3] shows the spectra and Figure |4] shows the corre- 
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Fig. 5. — Histograms of projected physical separations on the 
sky measured between double [O III] A5007 emission components 
in the long-slit spectra for the entire sample (black solid), objects 
with spatially compact emission components (red dashed), and ob- 
jects with spatially extended emission components (blue dotted). 
The median physical separation is 1.1 h^ kpc. See Section l4~2l and 
Figures [T] and [2]for a description and examples of spatially compact 
and spatially extended emission. The distribution of spatial sep- 
arations indicates that the double-peaked AGN emission lines are 
produced by ~kpc-scale dual AGNs and/or ~kpc-scale outflows, 
jets, or rotating gaseous disks 

sponding SDSS images of these galaxies. 

SDSS J 080740. 99+390015. 2 

Althou gh this galaxy was classified as a double-peaked 
AGN in IWang et all ([2009D . its spectra display triple- 
peaked narrow lines and single-peaked broad lines (Fig- 
ure [3j top), a line profile that could be explained by 
outflows or jets that affect the NLR only. We exclude 
this galaxy from the analysis in this paper, as we desire 
a homogeneous sample of double-peaked AGNs. 

SDSS J 00091 1 . 58- 003654 • ? 

For observations of this galaxy, we purposefully aligned 
one of the slit position angles (9 = 67°. 0) so that it 
would include the nearby spiral galaxy to the southwest 
of the target (Figure 2J middle). This spiral galaxy has 
photometric redshifts from SDSS of £ p hoto = 

0.083 ± 

0.017 (measured with the template fitting method) and 
%>hoto = 0.066 ± 0.022 (measured with a Neural Net- 
work method), which are similar to the target galaxy's 
spectroscopic redshift z = 0.073. If the spiral galaxy is 
indeed at a similar redshift as the target galaxy, then the 
spatially extended emission we observe (Figure [3j mid- 
dle) could be H/3. We measured the wavelength of the 
centroid of the emission (as in Section l4TT]h and if the 
emission is indeed H/3 then this wavelength corresponds 
to a spectroscopic redshift z = 0.071 or a velocity offset 
600 km s _1 blueward of the target galaxy. This is sug- 
gestive of the type of rich environment that we expect 



for dual AGNs, which are created from galaxy mergers. 

SDSS Jl 00145. 30+283330. 3 

This galaxy has an interacting companion to the north- 
west (Figure IU bottom), and we aligned one of the slit 
position angles (9 = 142°. 3) to include this companion. 
The long-slit spectra show emission features physically 
linking the two galaxies (Figure [3j bottom), possibly in- 
dicating an inflow or outflow of gas between the two 
galaxies. This could be an indication of the triggering 
mechanism that funnels gas onto SMBHs during galaxy 
mergers. 

7. RESULTS AND DISCUSSION 

7.1. All Double-peaked AGNs We Observed Have 
Double Emission Components with 
Kiloparsec-scale Separations 

We find that all of the double-peaked AGNs in our 
sample exhibit two spatially distinct AGN emission com- 
ponents, with projected separations on the sky ranging 
from hundreds of pc to a few kpc (the full range is 0.2 
kpc < Ax < 5.5 hjQ kpc). The median separation 
is 1.1 hjQ kpc, and the distributions of separations are 
similar for the subsamples with spatially compact and 
spatially extended emission components (Figure [5]). We 
note that, given the 3" diameter of an SDSS fiber, emis- 
sion components with separations > 3" (which, for our 
highest-redshift object, corresponds to Ax > 15 K^q kpc) 
will not fall within a single SDSS fiber. 

The range of projected separations we measure indi- 
cates that the mechanism or mechanisms producing the 
double-peaked AGN emission lines are neither a very 
small scale nor a very large scale effect. Hence, the 
data are not explained by the rotation of gas in disks 
of sizes <100 pc, bulk flows of gas on small scales near 
the galaxy's center, small-scale outflows, or AGN pairs 
or AGN outflows with ~10 kpc separations. Instead, the 
spatial separations we measure indicate that the double- 
peaked AGN emission lines in our sample are produced 
by ~kpc-scale dual AGNs and/or ~kpc-scale outflows, 
jets, or rotating gaseo us disks. 

We note that ISmith et all (|2Q12j ) identify 30 
double-peaked AGNs with peaks that have symmet- 
ric profiles or are nearly equal in flux and sug- 
gest that they are the products of rotating gaseous 
disks. Six of our objects are included in their 
sample of 30: SDSS J081542.53+063522.9, SDSS 
J084049.47+272704.8, SDSS J124813.82+362423.6, 
SDSS J130724.08+460400.9, SDSS J171544.05+600835.7 
(whic h X-ray observations in dicate is a dual AGN sys- 
tem; IComerford et all 120111 ). and SDSS J210449.13- 
000919.1. 

7.2. Physical Separation - Velocity Relation 

The correlations between the projected physical sep- 
arations and line-of-sight velocity separations of the 
double-peaked AGNs could help distinguish between dif- 
ferent sources of the emission, since dual AGNs may ex- 
hibit different correlations than NLR kinematics in single 
AGN. 

If a subset of our sample are dual AGNs, then the 
spatial separations and velocity offsets we measure are 
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Fig. 6. — Relation between projected physical separations on 
the sky measured between double [O III] A5007 emission compo- 
nents in the long-slit spectra and line-of-sight velocity separations 
between double-peaked [O III] A5007 emission lines in the SDSS 
spectra. The objects with spatially extended emission in the long- 
slit spectra (open squares) have a Spearman's rank correlation co- 
efficient of 0.34 (where the significance level of its deviation from 
zero is 0.05), suggesting they may be preferentially produced by 
AGN outflows. In contrast, the objects with compact emission in 
the long-slit spectra (filled circles) have no significant correlation, 
suggesting they may be preferentially produced by dual AGNs. 

snapshots of the dual AGN dynamics in the stages af- 
ter the initial galaxy merger but before SMBH coales- 
cence. For a fixed, three-dimensional physical separation 
between dual AGNs, we might naively expect a larger 
(smaller) line-of-sight velocity for a larger (smaller) pro- 
jected physical separation. When the orbiting dual 
AGNs are oriented almost along the line of sight, the 
projected physical separation will be small and the bulk 
of the velocity will be transverse, resulting in a small 
line-of-sight velocity measurement. In comparison, dual 
AGNs with a range of three-dimensional physical sep- 
arations may wash out any relation between projected 
physical separation and line-of-sight velocity. 

A positive correlation between velocity and physical 
separation could also result from AGN outflows, where 
a higher- (lower-) velocity outflow enables ionized gas 
to emit [O III] out to larger (smaller) distances from 
the central AGN. In contrast, the opposite correlation 
should occur for AGN jets. If a jet is oriented almost 
along (transverse to) the line of sight, we should measure 
a large (small) line-of-sight velocity but a small (large) 
projected physical separation. 

For the full sample, we plot the projected physical sep- 
arations between the double [O III] A5007 emission com- 
ponents seen in the long-slit spectra against the line-of- 
sight velocity separations between peaks in the double- 
peaked [O III] A5007 emission lines in the SDSS spectra, 
and we find that the Spearman's rank correlation coef- 
ficient is 0.16 and the significance level of its deviation 



Fig. 7. — Relation between projected physical separations on the 
sky measured between double [O III] A5007 emission components 
in the long-slit spectra and the difference AO between the measured 
position angle of the emission on the sky (from long-slit spectra) 
and the isophotal position angle of the major axis of the object 
(from SDSS r-band photometry). Since AO is only a meaningful 
measurement for galaxies with some degree of ellipticity, we plot 
only the 44 of our 81 galaxies that have ellipticities e > 0.1. Of 
these, we find that 37^g 2 % of the objects with spatially extended 
emission in the long-slit spectra (open squares) and 28^° % of 
the objects with spatially compact emission in the long-slit spectra 
(filled circles) have AO values that are, within 3cr, consistent with 
zero. These fractions are more than twice the expectation for a 
uniform distribution of AO between 0° and 90°, indicating that 
the double AGN emission components are preferentially produced 
in the plane of the host galaxy. 



from zero is 0.15 (Figure [6]). This significance level means 
that there is a 15% probability of observing a correlation 
this strong or stronger in uncorrected data purely by 
chance. 

However, the correlation strengthens to a Spearman's 
rank correlation coefficient of 0.34 (significance 0.05) for 
the subsample with extended spatial emission and essen- 
tially disappears for the subsample with spatially com- 
pact emission (correlation coefficient -0.04, significance 
0.77; see Section 14.21 for a discussion of extended and 
compact emission). The physical separation - velocity 
relation thus seems to separate into two populations: the 
spatially extended AGN systems and the spatially com- 
pact AGN systems. 

The spatially extended AGN population's correlation 
of increasing projected physical separation with increas- 
ing line-of-sight velocity is inconsistent with our expec- 
tations for AGN jets. However, the correlation might be 
explained if the population's spatially extended emission 
features are preferentially produced by AGN outflows. 
The correlation also matches that for dual AGNs at fixed, 
three-dimensional separations, but NLR effects would be 
necessary to explain the spatially extended emission. 

The spatially compact AGN population's lack of cor- 
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relation between projected physical separation and line- 
of-sight velocity could be explained by dual AGNs with 
a variety of three-dimensional physical separations. The 
spatially compact emission is consistent with the scenario 
of two compact AGNs with no extended NLR features. 

Although the physical separation - velocity relations 
hint that the spatially extended AGNs may be prefer- 
entially produced by outflowing AGNs and the spatially 
compact AGNs may be preferentially produced by dual 
AGNs with a range of three-dimensional physical sep- 
arations, we stress that additional observations are re- 
quired for confirmation of dual AGNs and outflowing 
AGNs. Further, the subpopulations of double-peaked 
AGNs with spatially compact or spatially extended emis- 
sion are likely not clearly defined as the products of either 
dual AGNs or NLR kinematics, as the emission profile of 
a double-peaked AGN could be produced by a combina- 
tion of both effects. 

7.3. Double Emission Components Are Preferentially 
Aligned with the Host Galaxy Major Axis 

Our measurements of the position angles of the dou- 
ble emission components can be used to determine their 
orientations in the host galaxies. We define AO as the 
difference between the isophotal position angle of the ma- 
jor axis of the galaxy (which is determined from SDSS 
r-band photometry) and the position angle of the two 
AGN emission components on the sky (which we mea- 
sure as # s ky in Section I47T]) . However, the isophotal po- 
sition angle of the major axis of the galaxy, and hence 
A#, is only a meaningful measurement in systems where 
the photometry shows the galaxy has some degree of el- 
lipticity. Consequently, we focus only on the 44 of the 
81 galaxies in our sample that have ellipticities e > 0.1 
(where ellipticities are computed in Section CT2]) . 

For these 44 galaxies, we find that the position an- 
gle of the two AGN emission components on the sky is 
often aligned with the isophotal position angle of the 
major axis of the object, such that the difference AO is 
close to zero (Figure Ej). In fact, 32±| % (14 objects) of 
the full sample of 44 galaxies, 37^g 2 % (7 objects) of the 
subset with spatially extended emission, and 28±\° % (7 
objects) of the subset with spatially compact emission 
have AO values that are within a 3a deviation from zero. 
These fractions are more than twice the expectation for a 
uniform distribution of AO: if the errors on AO remained 
the same but the AO values were distributed uniformly 
between 0° and 90°, only 14^% of the sample would 
have AO values within a 3a deviation from zero. 

Expressed another way, 30+g % of the 44 galaxies have 
AO < 10°, which is almost three times the expectation of 
lltl % for a flat distribution of AO between 0° and 90°. 

In addition, 55 tg % of the 44 galaxies have AO < 20°, 
which is more than twice the expectation of 22^ % for 
a flat distribution of AO between 0° and 90°. 

Alignment of the two AGN emission components with 
the major axis of the object (or low AO) indicates that 
the emission is produced in the plane of the galaxy, 
which is consistent with emission produced by dual 
AGNs (whose orbits are driven by the potential of the 
host galaxy and dyna mical friction from the host galaxy 
stars; in simulations, iBlecha et al.ll2Q12l also see align- 
ment of the dual AGNs with the ellipticities of the host 



gala xies) and some AGN-driven outflows and je ts (see, 
e.g. JWhittle fc Wilsoiill2QQ3 iRosario et al.ll2010h . Other 
AGN-driven outflows and jets can be more orthogonal 
to the host galaxy plane, appr oaching AO ~ 90° (e.g., 
IWhittle et aLlll988t iCecil et al.lkOOll ). but we do not an- 
ticipate many dual AGNs at these orientations. Based 
on these expectations, low AO could be used as a selector 
for double-peaked AGNs that are, with additional follow- 
up observations, more likely to yield detections of dual 
AGNs than double-peaked AGNs with high AO. 

We also plot AO against the line-of-sight velocity sepa- 
rations between [O III] A5007 emission components, but 
we find no significant correlation between them. 

7.4. The Most Promising Dual AGN Candidates 

While all double-peaked AGNs can be considered can- 
didate dual AGNs, we attempt to identify observational 
signatures that select the most promising candidates. We 
find that two spatially compact AGN emission compo- 
nents, both emission components in the plane of the host 
galaxy, and two stellar nuclei spatially coincident with 
the two emission components may be indicators of dual 
AGNs. In this way, we can focus limited observational 
resources for follow-up on the double-peaked AGNs that 
are more likely to yield dual AGN detections. We report 
these promising candidates in Sections 17.4.11 and 17.4. 2\ 
but definitive proof of dual AGNs requires the spatial res- 
olution of X-ray or radio cores associated with each AGN 
(for details see IComerford et all l201ll iFu et al.l [20111 
iTingav Waythll2011h . 

7.4.1. Fourteen Double-peaked AGNs with Emission in the 

Plane of the Host Galaxy 

Since the dynamics of kpc-scale separation dual AGNs 
are dominated by the host galaxy potential and dynam- 
ical friction from the host galaxy's stars, dual AGNs 
should preferentially reside in the plane of the host 
galaxy. In Section \7\3\ we used this expectation to iden- 
tify 14 promising dual AGN candidates that have double 
optical AGN emission components oriented in the plane 
of the host galaxy, where the difference in position angles 
AO are within 3a deviations from zero. These 14 candi- 
dates include seven double-peaked AGNs with spatially 
extended emission (SDSS J084049.47+272704.8, SDSS 
J095207.62+255257.2, SDSS J160524.59+152233.5, 
SDSS J163316.03+262716.3, SDSS J210449.13-000919.1, 
SDSS J225510.12-081234.4, SDSS J230442.82-093345.3) 
and seven double-peaked AGNs with spatially com- 
pact emission (SDSS J011659.59-102539.1, SDSS 
J080418.23+305157.2, SDSS J093024.84+343057.3, 
SDSS J102325.57+324348.4, SDSS 

J144804.17+182537.9, SDSS J162939.58+240856.0, 
SDSS J225420.99-005134.1). The seven double-peaked 
AGNs with spatially compact emission may be the most 
promising dual AGN candidates of the 14, since we 
showed in Section 17.21 that spatially compact systems 
may be preferentially produced by dual AGNs. 

7.4.2. Three Double-peaked AGNs with Double Stellar and 

AGN Emission Components that Are Potentially 
Coincident 

Of the seven double-peaked AGN in our sam- 
ple that exhibit double stellar components in AO 
imaging (Section [5]), we find that three (SDSS 
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J095207.62+255257.2, SDSS J115523.74+150756.9, and 
SDSS J123915.40+531414.6) have double AGN emission 
components with spatial separations and position angles 
comparable to the spatial separations and position an- 
gles between the double stellar components in the corre- 
sponding AO image (Figured]). In each case, the spatial 
separations are consistent to within our la errors and the 
position angles are consistent to ~ 15°, suggesting that 
the double AGN emission components are likely coinci- 
dent with the double stellar components. This scenario 
is a signature of a late-stage galaxy merger in which the 
stellar remnants of the original progenitor galaxies ac- 
company the two AGNs as they evolve into a kpc-scale 
pair. 

Consequently, these galaxies are three of the most com- 
pelling dual AGN candidates in our sample. One of these, 
SDSS J095207.62+2 55257.2, was al s o iden tified as a dual 
AGN candidate by iMcGurk et all (|2011ft and iFu et al.l 
(2012), who both deem it a confirmed dual AGN sys- 
tem based on integral-field spectroscopy that shows AGN 
emission coincident with double stellar components in 
the AO imaging. While double stellar nuclei coincident 
with double optical or near-infrared AGN emission com- 
ponents is strong evidence for dual AGNs, such a config- 
uration could also be explained by a single AGN between 
the two continuum sources. In addition, dual AGNs 
may not always have double stellar components visible 
in AO, as is the case for two systems that are strong 
dual AGN candidates based on X-ray and other obser- 
vations but exhibit onl y a single source in A O imaging 
(|Comerford et al.l 120111 : IB arrows et al.l l2012f ) . In cases 
like these, the second stellar component could be hidden 
by dust obscuration or could be too faint to be within 
the detection threshold. While the three galaxies dis- 
cussed here are good candidates for dual AGNs, addi- 
tional follow-up observations are required to confirm or 
refute them as dual AGNs. 

We note that SDSS J123915.40+531414.6 is an unusual 
case because the two emission features corresponding to 
the double peaks in the SDSS spectrum have a signifi- 
cantly smaller spatial offset (1.3 h^Q kpc, or / /39) than 
the spatial offset between the two stellar compo nents vis- 
ible in the AO image (4.1 kpc, or 1//24; IFu et al.l 
I2011al) . Instead, we find a third emission component (la- 
beled B in Figure [8j bottom) between the blueward (la- 
beled A) and redward (labeled C) features corresponding 
to the double peaks. This third emission component has 
a spatial offset (4.17 ± 0.72 kpc, or l'/20 ± 0'/21) 
and position angle (32°. 5 ± 3°. 3) relative to the blue- 
ward component that is consistent with those of the dou- 
ble stellar components. Note that there is also a fourth 
emission component (labeled D), visible only in the two- 
dimensional spectrum taken at position angle = 124°. 4, 
that is ~ 5" from the target galaxy's stellar continuum 
and may be reflective of an outflow. The long-slit and 
AO observations of SDSS J123915.40+531414.6 suggest 
it might host a triple AGN system or a dual AGN system 
with outflows. 

8. CONCLUSIONS 

We have obtained follow-up long-slit spectroscopy of 
81 double-peaked narrow-line AGNs in SDSS at 0.03 < 
z < 0.36 with Lick, Palomar, and MMT Observatories. 
The purpose of these observations was to obtain spatial 
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Fig. 8. — Same as Figure [TJ but for the three double-peaked 
AGNs that are some of the best dual AGN candidates. They 
are promising candidates because the spatial offsets and position 
angles of their emission features are consistent with those of the 
double stellar components visible in AO imaging. The A, B, C, 
and D labels in the bottom long-slit spectra show the positions of 
the four emission components found in SDSS J123915. 40+531414. 6 
(see Section [7. 4. 21 for details). 

information about the sources of the AGN emission, to 
aid in distinguishing whether the double-peak morpholo- 
gies are a result of dual AGNs, rotating gas in a disk, or 
AGN jets and outflows and on what physical size scales. 
We measure the projected spatial separation and posi- 
tion angle of the double AGN emission components on 
the sky, and compare to the orientation of the host galaxy 
in SDSS photometry as well as multiwavelength observa- 
tions of the double-peaked AGNs. Our main results are 
summarized below. 

1. In all of the 81 double-peaked AGNs we observe, 
the long-slit spectra reveal double AGN emission com- 
ponents with projected separations on the sky of 0.2 
kpc < Ax < 5.5 hjQ kpc (median Ax = 1.1 h^Q kpc; 
Figure [5]). These results indicate that the double-peaked 
AGN emission lines in our sample are produced neither 
by very small-scale gas rotation and outflows nor by ~10 
kpc scale AGN pairs or outflows. Instead, our observa- 
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tions are consistent with the double AGN emission com- 
ponents originating in ~kpc-scale dual AGNs or ~kpc- 
scale outflows, jets, or rotating gaseous disks. 

2. The double-peaked AGNs can be divided into two 
populations: those with double spatially extended emis- 
sion components whose line-of-sight velocity differences 
increase with projected physical separation on the sky 
(34/81, or 42^5% of our sample) and those with dou- 
ble spatially compact emission components whose line- 
of-sight velocity differences are independent of projected 
physical separation (47/81, or 58^ % of our sample). 
The physical separation - velocity trends shown in Fig- 
ure [6] suggest that the spatially extended sample may be 
preferentially produced by AGN outflows (since higher- 
velocity outflows can produce [O III] emission at larger 
distances) and the spatially compact sample may be pref- 
erentially produced by dual AGNs with a range of actual 
three-dimensional separations. Consequently, we suggest 
that double-peaked AGNs with double, spatially com- 
pact emission components may be fruitful targets for 
dual AGN searches. 

3. In 32^g % of the objects in our sample, we find 
that the orientations of the double AGN emission com- 
ponents on the sky align with the major axis of the host 
galaxy, such that the differences AO are within a 3<r de- 
viation from zero (Figure [7j). This result is more than 
double the expectation (14+4 %) for a uniform distribu- 
tion of AO. Since dual AGNs should orbit in the plane 
of the galaxy, small AO could be a useful indicator for 
promising dual AGN candidates. Our proposal of small 
AO as an indicator for dual AGNs is supported by sim- 
ulations that show a strong alignment between the ori- 
ent ations_oXdja^J_^AGNsand the major axis of the host 
galaxy (jBlecha et al.ll2012l ). 

4. We identify 17 particularly promising dual AGN 
candidates for follow-up observations. Fourteen double- 
peaked AGNs (listed in Section 17.4. lj) have double AGN 
emission components in the long-slit spectra that are ori- 
ented in the plane of their host galaxies, as expected for 
dual AGNs. In an additional three double-peaked AGNs 
(listed in Section l7.4.2l and shown in Figure[8]), the spatial 
offset and position angle we measure between the double 
AGN emission components are consistent with the spa- 
tial offset and position angle between the double stellar 
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components in AO imaging. This suggests that the dou- 
ble AGN emission components may be coincident with 
the double stellar components, which is a scenario ex- 
pected for dual AGNs. 

While the one-dimensional spectra of double-peaked 
AGNs alone do not reveal much about the sources of 
the emission, we have shown that the spatial informa- 
tion about the emission provided by follow-up long-slit 
observations can help distinguish between between gas 
rotation, AGN outflows and jets, and dual AGNs on dif- 
ferent physical scales. In this way, long-slit observations 
enable identification of the most compelling dual AGN 
candidates for confirmation via X-ray or radio observa- 
tions that spatially resolve two AGNs. 
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